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Abstract 
Modeling of biomass pyrolysis kinetics is an essential step towards reactors design for energy production. Determination of the 
activation energy, frequency factor, and order of the reaction is necessary for the design procedure. Coats and Redfern's work 
using the TGA data to estimate these parameters was the cornerstone for modeling. There are two significant problems with 
biomass modeling, the first is the determination of the kinetic triplet (Activation energy, Frequency factor, and the order of 
reaction), and the second is the quantitative analysis of products distribution. Methods used in modeling are either One-step or 
Multistep methods. The one-step techniques allow the determination of kinetic triplet but fail to predict the product distribution, 
whereas multistep processes indicate the product's distribution but challenging to estimate the parameters. Kissinger, Coats, 
and Redfern, KAS, FWO, Friedman are one-step methods that have been used to estimate the kinetic parameters. In this work, 
after testing more than 500 data points accessed from different literature sources for coal, oil shale, solid materials, and biomass 
pyrolysis using one-step global method, it was found that the activation energy generated by KAS or FWO methods are related 
as in the following equations: Egas = 0.9629 * Epwo + 8.85, with R? =0.9945 or Epwo = 1.0328 * Egas — 8.0969 with R?= 0.9945. 
The multistep kinetic models employed the Distributed Activation Energy Model (DAEM) using Gaussian distribution, which 
suffers from symmetry, other distributions such as Weibull, and logistic has been used. These multistep kinetic models account 
for parallel/series and complex, primary and secondary biomass reactions by force-fitting the activation energy values. The 
frequency factor is assumed constant for the whole range of activation energy. Network models have been used to account for 
heat and mass transfer (diffusional effects), where the one-step and multistep could not account for these limitations. Three 
network models are available, the Bio-CPD (Chemical Percolation Devolatilization) model, Bio-FLASHCHAIN, and the Bio-FG- 
DVC (Functional Group Depolymerization Vaporization Crosslinking models). These models tried to predict the product 
distributions of the biomass pyrolysis process. 
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Introduction 

Thermal decomposition reactions of biomass play a crucial role in its utilization processes (Balan, 2014). Biomass, which covers from 
tiny grass to massive trees and from animal waste to small insets (Basu, 2010), is a renewable energy source that produces heat and 
energy whenever needed (Anca-Couce, 2016; McKendry, 2002). Biofuels contain various pyrolyzing species with different reactivity 
and are probably influenced by surrounding species (Carpenter et al., 2014). Biomass, by enlarge, comprises cellulose, which 
constitutes 38 — 50% of the total mass, hemicellulose, 20-32%, lignin 15-25%, and 5-13% other inorganic and extractive species 
(Hameed et al., 2019). There are three principal routes to release heat and energy from biomass: Pyrolysis, direct combustion, and 
gasification (Cuoci et al., 2007). The transformation of biomass can be performed either by biochemical or thermochemical methods. 
In biochemical conversion, sugar platform biorefinery includes physicochemical pretreatment of the biomass, the enzymatic hydrolysis 
of the carbohydrates to a fermentable sugar stream, and finally, the fermentation of the sugars by suitable microorganisms to the desired 
molecules (Simone and Michael, 2015). 
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The main thermochemical processes are gasification, pyrolysis, combustion, hydrothermal liquefaction, and hydrothermal 
carbonization (Behrendt et al., 2008; Funke et al., 2010; Basu, 2010; Sharma et al., 2015). Processes that convert biomass into liquid 
fuels and other products proceed via the pyrolysis route. 

The pyrolyzed bio-oil might undergoes upgrading (Zhang et al., 2007; Demirbas, 2003a, 2003b) to improve properties and meet the 
end-user specifications. The kinetics of thermochemical pyrolysis is the essence of reaction rate, product distribution, and reactor 
design. Biomass pyrolysis is the process that is positively affected by a variety of parameters such as heating rate (Mishra and Mohanty, 
2018), temperature (Angin Dilek, 2013), gas flow rate (Tripathi, et al., 2016), biomass composition, residence time, moisture content, 
particle size (Koufopanos et al., 1989(a) 1991(b); Sheth and Babu, 2006; Ranzi, et al., 2008; Vinu and Broadbelt, 2012), and type of 
reactor (Lédé, 2012). The importance of biomass pyrolysis kinetics is emanated from the need for gasifiers and combustors reactors 
design (Patra and Sheth, 2015). Design and optimization of reactors require estimating the relevant parameters such as activation 
energy, frequency factor, and reaction order that pertain to the decomposition process (Huang ef al., 2013; Meng et al., 2013). 
Numerous factors affect the pyrolysis rate, the yields, composition, and properties of the product classes. Temperature, pressure and 
heating rate, heating time (Gronli and Melaaen, 2000) affect the pyrolysis product distribution and are the main operating parameters 
(Newalkar et al., 2014). In addition, biomass properties (chemical composition, ash content, composition, particle size and shape, 
density, moisture content (Demirbas, 2004; Mehrabian, et al., 2012), heating value, etc.) also play an essential role. 

Pyrolysis involves a breakdown of large complex molecules of the biomass into several smaller molecules. The cellulose, 
hemicellulose, and lignin undergo pyrolysis differently, contributing to yields. Hemicellulose, which is a short branched high molecular 
weight linear polymer of B-(1—4)-D-glucopyranose units of different monosaccharides linked to each other by (1-4)-glycosidic bonds 
(Mohan et al., 2006) undergoes exothermic pyrolysis at temperatures between 150—350°C then, shifts to endothermic pyrolysis at 
higher temperatures (Broido and Nelson, 1975). Pyrolysis of hemicellulose yields more noncondensable gases such as CO2, CO, Hz, 
CH4, C2H2, C2H4, C2H6. In addition to bio-oil and char. Cellulose is made of unbranched fibrils composed exclusively of glucose. 
Thermal destruction of cellulose during pyrolysis proceeds through gradual degradation, decomposition, and char formation at slow 
heating and low reaction temperatures. Rapid volatilization is accompanied by levoglucosan formation at high pyrolysis temperature 
and fast heating rates. The endothermic cellulosic degradation reactions occur in the range of 275-350°C; while exothermic reactions 
take place at higher temperatures. Cellulosic reactions include hydrolysis, depolymerization, oxidation, decarboxylation, and 
dehydration. The main liquid products are levoglucosan-pyranose, glycolaldehyde, formic acid, anhydrous xylopyranose, 5- 
hydroxymethyl furfural, and 2-furaldehyde, etc. (Patwardhan et al., 2009) in addition to char and gaseous products. Lignin is a high 
molecular weight amorphous, a phenolic compound with different alcohol building units such as sinapyl, coniferyl, and p-coumaryl. 
Pyrolysis of lignin proceeds mainly in the temperature range between 200 and 500°C (Hosoya et al., 2008). The pyrolysis of lignin 
culminated in the production of unsaturated side-chains, including alcohol, aldehyde, guaiacol, stilbene, gases, char. 

There are a numerous studies on biomass pyrolysis kinetics for developing various kinetic models (Hameed et al., 2019; Lam et al., 
2012; Hu et al., 2007). Most pyrolysis kinetic models of biomass are considered lumped models since the kinetics are based on the 
yield of lumped products (i.e., gas, bio-oil, and solid-char). More light will be shed on the different modeling kinetics and review 
mathematical development in this research work. 


1 Materials and Methods 
1.1 Chemical Aspects 
1.1.1 Mechanism and reaction modeling 


Each of the structural constituents of biomass (cellulose, hemicelluloses, lignin, and extractives) is known to pyrolyze at different rates 
and by different mechanisms and pathways at various temperatures. These reactions are spread widely over the 150-550°C temperature 
range. The rates of reactions and extent of degradation of each of these reacting components depend on the process parameters, 
temperature, biomass type, particle size, heating input rate, and total pressure of the system. Bradbury et al., 1979, studied cellulose 
pyrolysis and reported series (incubation period to convert inert materials to active components) and parallel reactions in which the 
first stage reactions were characterized by high activation energy. The pyrolysis products are gases of similar elements with a different 
distribution, liquids of varying structure and types, and char solids of other formats. (Chen, et al., 2019; Yang et al., 2007; Kawamoto, 
2017). Chemical kinetics play a vital role in elucidating the characteristics of pyrolysis reactions and developing mathematical models 
needed for reactor design. 

The early use of thermogravimetric data (TG) to evaluate kinetic parameters of solid-state reactions was used by Coats and Redfern in 
1964. These authors advocated using a small size sample to eliminate the heat effect due to endo- or exothermic heat of reactions and 
maintain a linear heating rate and accurate temperature measurement. 

It is stated that for solid-state chemical reactions, the rate of decomposition of reactant is proportional to its unreacted concentration 
and reaction temperature, as described in equation (1): 
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<= kf (a) = k(1 — a)” (1) 


where a is the fraction of solid materials decomposed defined as indicated in equation (2) at any t , n is the reaction order, 


a= (wo —w,)/(Wo — Wr) (2) 
where w, is the initial weight of biomass, w, is the weight of biomass at any run time, t and wẹ is the final weight of biomass at the 


end of the run, k is the rate constant given by Arrhenius equation (3): 


k = k,exp(— E/RT) (3) 
in which k, is the frequency factor (s~*), E is the activation energy, (kJ /kmol), R is the universal gas constant, (8.314 kJ//kmol K) 
and T is reaction temperature in (K). 
Introducing the linear heating rate, 8 = dT /dt into equation (1) to obtain equation (4) 


da/dT = (k,/B) exp(—E/RT) (1 — a)” (4) 
If equation (4) is integrated, we obtain equation (5) 


a da Ko rT E 
Se = BS exp (-) ar (5) 


The RHS of equation (5) has no exact integral, but with some mathematical manipulations by substituting u = E/RT and using the 
relation of equation (6) (Coats and Redfern, 1964) to obtain equation (7) below: 


[exp ud = u!Merp™* Fp IO © 
Now, equation (5) becomes 
1-0-4)" ain = a = = E 
= —[1-Flex-z) 7) 
If the logarithm of Bai sides i is iken we obtain equation (8) as follows 
1-(1-a)1 =") _ KoR _2RT] E 
Inf (1-n)T? } =In (=) [1 E | RT (8) 


If reaction order in equation (5) is taken to be unity (a special case of equation (8)) and after taking the logarithm, equation (9) is 
obtained: 


-In(1-a@) 2RT 
in|] = ng a (9) 
Plotting the left hand Pe of ae (9) against the inverse of reaction temperature (1/T) would results in activation energy 


calculations from the slope of the straight line resulted from plotting the data. Simultaneously, the frequency factor could be determined 


from the intercept of the line with the y-axis. The quantity m*a - = is assumed constant (Coats and Redfern, 1964). 


Furthermore, the magnitude of the quantity 2RT/E within the range of biomass reactions, temperatures are relatively small compared 
with unity. Hence, it could be neglected; as a result of both cases, the intercept would value the frequency factor. Equation (9) has 
been extended or modified and developed further to model solid-state decomposition reactions. In the present work, several kinetic 
models based on equation (9) will be discussed and used by investigators within the scope of biomass pyrolysis kinetic modeling. The 
majority of kinetic models used by researchers employed a one-step global kinetic model or more than one-step kinetic models such 
as parallel and competitive reactions. The detailed lumped kinetic model, models with secondary bio-oil cracking, distributed activation 
energy model (DAEM) based on Gaussian distribution (Braun and Burnham, 1987), Weibull (Lakshmanan and White, 1994) and 
Gamma distributions (William and Tarutis, 1994), and nucleation growth model. Figure 1 shows the various kinetic models employed 
by investigators. As indicated, one-step and multistep models have been used. 
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Kinetic Modeling 


One-step global kinetic model 
(Primary reactions, No 
Mechanism 


One-step global kinetic models: Parallel and Competitive 
represent the conversion of reactions Models 


biomass to the volatiles and char 
as a first order single-step 
reaction. Kissenger, KAS, FWO 


and other Models are used to : 
estimate kinetic parameters or Avrami-Erofeev or 
experimentally determined. Prout-Tompkins 


Fig. 1 Various kinetic models employed for biomass reactions. 
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1.2 Kinetic modeling 
1.2.1 One-Step global kinetic model 


The one-step global kinetic models represented by a first-order single-step reaction are the simplest pyrolysis models representing the 
conversion of biomass to the volatiles (vapors and non-condensable) solid-char. The pyrolysis products only consist of char and only 
volatiles (White et al., 2011). This one-step kinetic model was also employed for intrinsic kinetics and under the influence of heat and 
mass transfer limitations (Antal et al., 1998; Prakash and Karunanithi, 2008). These models do not account for the reaction mechanism 
concerned with immediate reactions for volatiles and solid-char as final products. Several models such as Kissinger Model, Flynn- 
Wall-Ozawa (FWO), Kissinger-Akahira Sunose (KAS), and Friedman, Coats, and Redfern have been developed and used (Papari and 
Hawboldt, 2015). 


1.2.2 Kissinger model equation 


The widely used Thermal Gravimetric Analysis/Differential Thermal Analysis (TGA/DTA) for conducting experimental work is a 
thermal analysis in which the reacted mass is measured as a function of time at non-isothermal conditions with a constant linear heating 
rate. When the DTA data are plotted against the pyrolysis temperature, the curve exhibits a maximum mass loss rate at a specific 
temperature. At a maximum rate of decomposition at the peak of DTA curve, the first derivative of equation (1) would be zero, and 
the corresponding temperature Tm is termed characteristic temperature. Taking the second derivative at constant linear heat rate B = 


dT /dt equation (10) is obtained as follows: 


ae $e ( 5 4 LO kep 5) (10a) 
(a) gt ae ker (a) = (10b) 


Since the rate of conversion, i.e. da/dt is a non-zero value, then the square bracket term must be zero leading to Kissinger’s equation 
(Kissinger, 1957) indicated as equation (11) in which after taking the logarithm and rearrangement becomes: 


in (4) =In [=< eed), 2. (11) 


E da RTm 
If the reaction rate order is assumed unity, and for a given value of the linear heating rate 8, equation (11) is employed to estimate the 


frequency factor and activation energy by plotting[/n (E) against the inverse of several values of Tm- 


1.2.3 Flynn-Wall-Ozawa (FWO) model equation 


The right hand side (RHS) integral of equation (5) can be expressed through (Doyle, 1961) tabulated values as indicated in equation 


(12): 
aio exp (—s) aT = 250 (5) (12) 
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If the quantity 5 is larger than 20, the quantity p (=). can be approximated and become exact (Doyle, 1962) as in equation (13) 


which is obtained from the tables of the exponential values. 


log(E/RT) = —2.315 — 0.4567 E/RT (13) 
Whereas equation (13) becomes inexact if the quantity (=). is less than 20. Furthermore, equation (5), after introducing the 


approximation, becomes equation (14): 


So “5 = o(@) = 2S) exp(-E) ar = 22 p (5) = Game (14a) 

Infg(a)] = In (=) — 2.315 — 0.4567 = (14b) 
Rearranging: 

In(B) = in(= Eko 2 2.315 — 0.4567 — (14c) 


Now, for any evolved biomass loss quantity, i.e. a given constant value of the conversion function, a, the left-hand side of equation 
(14b) is independent of heating rate. Suppose the conversion function a is kept constant for several heating rate values. In that case, 
the LHS plot of equation (14c) against the pyrolysis temperature inverse will culminate in activation energy and frequency factor 
determination. FWO method's advantage is that it does not require any form of kinetic model other than the Arrhenius equation for 
temperature dependence. The FWO method is a model-free that involves measuring the temperatures corresponding to fixed values of 
conversion function a from experiments at different heating rates, B. 


1.2.3 Kissinger-Akahira Sunose (KAS) 


The KAS method is based on equation (5) with the approximation in equation (6) to obtain equation (15) 


a da =" __ KoRT? 2RT 
Ss sa go ae) oT ge (15) 


and if the nate = is << 1, then equation (15) can be written as indicated by equation (16) i.e. 


a da E 
[i = g(a) = exp (- 5) (16a) 
Rearranged after taking ma Daaa of both sides: 
Rko E 
in(4 5) = =In n() i RT (16b) 


The activation energy and frequency factor are determined if the LHS of equation (16b) is plotted against the inverse of pyrolysis 
temperature for several values of the heating rate at a constant conversion function g (æ) value. The difference between KAS and FWO 


is the method of evaluating and approximating the value of the integral of ( f z exp (- =) dT). Figure 2 shows the direct relationship 


between the activation energy obtained through KAS method and the activation energy estimated using FWO method. The closeness 
of the values predicted by both ways is related to the methodology of approximating the RHS, i.e. the integral part of the equation 
(15). The presented data in Fig. 2 are sourced from different studies and arbitrarily selected (more than 28 references) investigation of 
biomass, lignite coal and oil shale pyrolysis studies (Aboulkas and El-Harfi, 2008; Zou et al., 2010; Damartzis et al., 2011; Al-Ayed, 
2011; Katarzyna et al., 2012; Mishra and Bhaskar, 2014; Kongkaew et al., 2015; Heydari et al., 2015; Islam et al., 2015; Maurya et 
al., 2016; Ahmad et al., 2017; Chandrasekaran et al., 2017; Mishra and Mohanty, 2018; Xu Bang et al., 2018). As it can be inferred 
from the plots, Figure 2 and Figure 3, the activation energies estimated by KAS and FWO methods are related as indicated by equation 
(16). Equation (16a) predicts the KAS estimated activation energy as a function of the FWO calculated activation energy. Equation 
(17) is generated from actual experimental data taken from the pyrolysis of different biomass, coal, oil shale, and solid material 
reactions. 


Exas = 0.9629 * Epwo + 8.85, R2=0.9945 (17a) 
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Fig. 2 Comparison of the activation energies obtained for KAS and FWO methods from different literature. 
(sources: Mishra, et al., 2018 and refs. from Aboulkas and El-Harfi, 2008 until Chen et al., 2015) 
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Fig. 3 Comparison of the activation energies obtained for KAS and FWO methods from different literature sources. 
(sources: Mishra, et al., 2018 and refs. from: Aboulkas and El-Harfi, (2008) until Chen et al., 2015) 


When the plotted data are exchanged on the axis, the FWO estimated activation energy is given as a function of KAS estimated 
activation energy as indicated in equation (17); it is worth noting that this equation is generated through trend line of the data plotted 


on Fig. 3 rather than a rearrangement of equation (17a). 


Erwo = 1.0328 * Egas — 8.0969, R?=0.9945 (17b) 

The coefficient of determination R? the value which measures the goodness of the fit, is relatively high, indicating an excellent 
correlation between the KAS and FWO methods of predicted activation energies. Coefficients of equation (17) could be improved 
furthermore if more published data are included in the plot to become universal, although more than 500 data points have been used 
to plot Fig. 2 and 3. It is quite possible to generate the KAS activation energy values from the corresponding FWO calculated energy 
values and vice-versa for any biomass, coal, oil shale, and any other solid material without conducting actual calculations. It can be 
concluded that KAS method is basically the same as FWO method and vice-versa since the small difference between the results 
inferred between both ways is related to the approximation procedure. This newly obtained equation (17) can replace the process of 
calculating the activation energy of either KAS or FWO calculations work if either of the pertaining values of KAS or FWO are 
calculated from the experimentally found data. Also, the generated equation could be used to create a check value on the experimentally 
determined activation energy. It should be noted that equations (17a) and (17b) are obtained directly from plotting data. 
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1.2.4 Friedman method 


The un-widely used power-law method of Friedman (Friedman, 1960) is a technique for obtaining kinetic parameters that describe the 
thermal degradation using TGA data. This is a differential isoconversional method where it is possible to determine the activation 
energy of specific processes without knowing the kinetic equation. If the logarithm of both sides of equation (1) is taken, equation (18) 
is obtained as follows: 


In (5) =- + In(kof(a)”) (18) 
In this method, the conversion functions f (æ)is assumed constant, indicating that biomass decomposition is independent of reaction 
temperature and depends only on the mass rate loss. Friedman believed that there was a single value of the reaction order (n), which 
indicates that the value of E/R is the same for all levels of conversion. A plot of the LHS of equation (18) against the inverse of 
pyrolysis temperature would lead to the activation energy. However, the slope's actual value is not the same for all conversion levels 
(Mishra and Mohanty, 2018). Different forms of the conversion function are assumed, as shown in table 1. 


1.2.5 Coats-Redfern method 


The Coats-Redfern equation is an integral model-free method (Coats and Redfern, 1964), broadly used to calculate reaction order and 
pre-exponential factor (Damartzis et al., 2011). Equation (7) can be rewritten as in equation (18) 


a da koRT? 2RT E 
g(a) = ff OR = ET i — a Cr) (19) 
Using an asymptotic approximation for the resolution of Eqn. (19) and if (2RT/E<1), equation (20) can be obtained: 
g(a)\ _ RKo M E 
In (=) =In (Fe) RT (20) 


If equation (20) is compared with equation (16b), we note both equations are the same except exchanging the location of the quantity 
g(a) with that of the heating rate, £ . It is concluded that the three methods, namely, KAS, FWO, and Coats and Redfern, are basically 
the same, and the only difference is the type of model used in the Coats and Redfern equation. In other words, KAS method is a 
particular case of the Coats and Redfern equation since several models (in addition to reactants nature and composition, pyrolysis 
products distribution influenced by the operating conditions during the reaction) can be used to find the g(a) as indicated in Table 1. 


Table 1 List of reaction models used by researchers. 


Reaction Model Symbol f(a) g(a) 

Power law 1 403 3 a! 
p7 
Power law n=} p3 3 a 2a73 
Power law n=> pz 2 a a7? 
3 3 2) 4 Ve E- 
Power law n=> pZ 52? ~ za? 
One-D diffusion (Parabola law) pt Zani -ar 
Mampel (n = 1 ) Ft 1-a -1 
Mampel (n‘" order ) F” (1 -—a)” =n(1- a)" 
Avrami-erofeev (n = 4 ) At 4(1 —a)(—InQ1 — a))# (4n(1 — a) + 3)(-In(1 — ayi 
Avrami-erofeev (n =3 ) AB 3(1 = a) (—ln(1 = a))3 —(3In(1 — a) F= a3 
Avrami-erofeev (n = 2 ) A? 2(1 — a)(-In(1 — a))2 Qna —a) + 1)(—In(@t — a))? 
3 3 3 1y\71 1 1 1 -2 

3-D (Jander model) D3 5 (1 —a)? (1 -~a- a3) z(a ES 2) (a me 1) G-a) 
Contracting sphere R? 31- ae -20-a)3 
Contracting cylinder R? 21- ae -~(—a)2 
2-D diffusion (Valensi model) D? (—In(1 — a))72 (-In(@l -a))?(a - 1) 
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1.2.6 Multistep kinetic models (Secondary reactions mechanism) 


The one-step global kinetic models represent the characteristics of the simple biomass conversion to gas, bio-oil, and char; it does not 
account for or consider the nature of the biomass, the constituents, the reaction operating conditions such as pressure, liquid hourly 
space velocity, rate of heat input, thermodynamics, etc. The simple one-step global kinetic models are used to estimate the kinetic 
parameters for immediate reactions but did not tackle the kinetics of secondary reactions that takes place in parallel/series or an 
elaborate manner after the initial decomposition of the biomass (Vinu and Broadbelt, 2012; Corbetta et al., 2013). The one-step global 
kinetic model's main advantages are the ease and direct determination of the activation energy, the frequency factor, and the needless 
of the reaction mechanism, and the disadvantage is the inability to predict the product's composition distribution. As a result, it is 
imperative to search for models that could predict product distribution, including primary and secondary biomass pyrolysis processes. 
Researchers formulated models to include the primary’ reactions products as secondary reactions to produces materials similar to the 
direct products, i.e., gas, tar, and char, through cracking and polymerization reactions (Thurner, 1981; Bradbury, et al., 1979; Scott, et 
al., 1985; Agrawal, 1988). Several kinetic schemes have been proposed for combining primary degradation of biomass and secondary 
decomposition of volatile products (higher molecular weight hydrocarbons). As early as 1975, Broido and Nelson formulated a 
multistep kinetic model by lumping the consecutive steps of cellulose reaction into a single, overall, first-order reaction as follows: 


Volatile tars 
Cellulose a 
Char + low molecular weight voaltiles 


Broido and Nelson's model was further reformulated (Bradbury et al., 1979) to include a cellulose activation reaction (incubation 
period). The cellulose is catalyzed through temperature treatment to initiate the reaction shown in Figure 4. The proposed scheme is 
presented as follows: 


Biomass (inactive cellulose) — Active cellulose, E = 242.672 kJ/kg 


Condensable volatiles 


Biomass (inactive cellulose) ———~ Active cellulose ait 


Char + gases (CO,CO,, H,O) 


Fig. 4 Bradbury and coworkers model of cellulose reaction. 


The model developed by Bradbury and coworkers (Bradbury et al., 1979) was further modified (Diebold, 1994) to consist of several 
intermediate and secondary reactions, as indicated in Figure 5. 


Secondary gas 


Cellulose Active Cellulose Primary vapors Secondary tar (v) 


Char (a) + H,O (a) 


Char (d)+ H,O(d) 


Fig. 5 Seven global reaction scheme for cellulose pyrolysis. 


Shafizadeh and Bradbury's model included both primary biomass pyrolysis reactions and secondary tar decomposition. In this model, 
(Bradbury et al., 1979), cellulose was first converted to an intermediate component (active cellulose), which is further decomposed to 
other secondary products. Furthermore, they assumed an initiation reaction for the cellulose to be functional (incubation period) in 
which no weight loss is observed. The apparent activation energy was estimated as 424.67kJ/kg. Two competitive first-order primary 
reactions to produce condensable volatiles and char cum gases (H20, CO, and CO2) were assumed to take place. The proposed model 
gave good predictions for the product yield (Di Blasi, 2008), but kinetic parameters’ determination was difficult. Researchers have 
used this model in its original form and extended formulation (Sheth and Babu, 2006; Varhegyi et al., 1994). 
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Kinetic models characterizing the main structures of devolatilization of the three components of the biomass have been developed by 
Cuoci and coworkers (Cuoci et al., 2007). Complete prediction of product and by-product has been reported based on a semi-detailed 
kinetic mechanism. This proposed mechanism included 15 reactions with 30 lumped product species. The development of Cuoci et 
al., was extended by Ranzi and associates (Ranzi et al., 2008) to a mechanistic model to describe the biomass pyrolysis, 
devolatilization, and gas-phase reaction of the product gas species from cellulose, hemicellulose, and lignin of the biomass. Other 
researchers (Corbetta et al., 2013) extended the works of Ranzi et al., to encompass gas-phase secondary reactions; further extension 
of this work was conducted (Couce et al., 2014) to include secondary char formation reactions under diffusional influences and in the 
absence of ash-alkali-metals effect to avoid sugar formation. It can be concluded from the published literature that all proposed 
mechanisms are the same if operating conditions such as temperature, pressure, and reaction environment, the residence time of 
products, biomass particle size, slow/fast heating rate, biomass constituents, and nature are alienated from influencing the biomass 
reactions. Figure 6 shows the different proposed mechanisms and the progress of suggested decomposition, depolymerization chemical 
reactions, and effect of reaction products. 


Volatiles (2) + Gases (2) 


Volatile (1)+ Gases (1) 


Gases + Operating conditions 


Char + Tar +Gases 


Fig. 6 Pyrolysis reactions and product variations (Varhegyi et al., 1997; Sheth and Babu, 2006; Ranzi et al., 2008). 


Biomass (active or inactive) 


Volatile (3) + Gases (3) 


During pyrolysis, the operating conditions, such as fast or low heating rate, could produce more or less tar than char (Klass, 1998). 
Fast heating yields higher volatiles and reactive char than the slow heating process (Basu, 2010). In contrast, slow heating with longer 
residence time produces more char due to secondary reaction between formed char and volatiles (De Groot and Shafizadeh, 1984). 
The Distribution of Activation Energy Model (DAEM) has been recently used to understand the pyrolysis process complexity. The 
DAEM considers many primary and secondary biomass independent reactions with the activation energy for each reaction estimated 
by the Gaussian distribution continuous function (Hu et al., 2016). Soria-Verdugo et al., (2014) found a difference between the 
activation energies and frequency factors calculated at low and high heating rates during the TGA tests. DAEM was used to model 
municipal solid waste and agricultural residues pyrolysis at low heating rates using non-isothermal TGA analysis (Bhavanam and 
Sastry, 2015). The DAEM could model and explain the devolatilization curve for microalgae pyrolysis (Ceylan and Kazan, 2015). A 
three pseudo-components DAEM was developed using probability density function (PDF) to describe the biomass conversion and 
thermal decomposition rates during Rhus Typhina pyrolysis (Liu et al., 2019). 

The Gaussian distribution disadvantage is that it is symmetric, whereas the real reactivity distributions are asymmetric (Burnham and 
Braun, 1999). The asymmetry can be accounted for using other distributions, e.g., the Weibull (Lakshmanan and White, 1994) or 
Logistic distributions (Cai et al., 201 1a, b). These models, multistep and one-step global kinetics, have been employing TGA tests and 
did not account for the mass and heat transfer diffusional limitations. At higher heating rates or large biomass particle size, the effect 
of mass and heat transfer become predominant. All the lumped product models (gases, tar, and char) were not able to explain the 
complex biomass pyrolysis reactions and the formation of several species. Therefore, network models have been utilized to address 
these diffusional limitations (Hameed et al., 2019). 


1.2.7 Network models 


Pyrolysis basically is a devolatilization process, which is a term for removing volatiles from solid materials such as coal, oil shale, and 
others. Devolatilization is similar to the desorption step in heterogeneous chemical reactions (catalytic reactions), i.e., one diffusional 
direction process (products desorption from inner of catalyst, reaction site to external surface of catalyst). Devolatilization generally 
involves reactions such as molecular-level depolymerization and polymerization (Wan et al., 2015). As shown in Figure 7, the 
diffusional process steps are divided into four stages in series and are considered the main stages of the network models (Borah et al., 
2011). 
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Stage 1 Stage 2 Stage 3 Stage 4 
r ; ? Products mass Diffusion through 
Solid material Intermediate - 
: : ohio transfer from pores to external 
Chemical reaction. formation, series 5 : 
: i reaction site to the surface of 
to give products. reactions : . 
surface material 


Fig. 7 Stages of solid material reaction and diffusion of products to external surface in network models. 


Researchers have developed three network models, which are also called structural models; the typical biomass network models are 
the Bio-CPD (Chemical Percolation Devolatilization) model (Sheng and Azevedo, 2002), the Bio-FLASHCHAIN model (Serio, 1997), 
and the Bio-FG-DVC (Functional Group Depolymerization Vaporization Crosslinking) model (Niksa, 2000). The bio-Chemical 
Percolation Devolatilization model (CPD) proposes that the solid materials are assumed to be clusters joined together through chemical 
bridges. Some of these connectors break at low temperatures, while others break at higher reaction temperatures. Low temperature 
breaking bridges give rise to gases and bio-oil (tar) formation, whereas those breaks at higher reaction temperatures give rise to char 
formation (Sheng and Azevedo, 2002). Workers (Lewis and Fletcher, 2013) predicted sawdust pyrolysis yields from a flat-flame burner 
using the CPD model. They assumed that the yield is the sum of the weighted average yield of the biomass's main components. Bio- 
FLASHCHAIN model was modified for biomass pyrolysis (Niksa, 2000) since initially it was developed for rapid coal pyrolysis. 
DISCHAIN, DISARAY, and FLASHTWO are the three concepts of the Bio-FLASHCHAIN theory. DISCHAIN concept is 
depolymerization, and cross-linking indicates the evolution of tar from the main reactant chains to small fragments; the produced 
monomers may react to the long pieces forming tar precursors by cross-linking to char, respectively. DISARAY is similar to 
DISCHAIN concept but applied to the two-dimensional structure. The FLASHTWO concept accounts for pressure effects, which has 
not been studied (Niksa, 1987). 

The third structural FG-DVC model consists of two implicit models, responsible for forming light gases due to the decomposition of 
the biomass functional groups and the polymerization of fragments as time proceeds to form tar. The first FG model accounts for the 
Function Group, and the second is the DVC (Depolymerization, Volatilization, and the Cross-linking model or the tar-formation model 
(Solomon et al., 1993). The DVC model determined the molecular weight distribution of char formed (Solomon and King, 1983). 


Conclusions 


The one-step global method is widely and most popular model used by researchers for determining the kinetic parameters. The lumped 
model is a free kinetic model. Multistep kinetic models use different reaction mechanisms and consider primary and secondary biomass 
reactions to assess the kinetic parameter cumbersomely. Distribution activation energy model using Gaussian and other distributions 
to estimate the kinetic parameters. Since none of the one-step and multistep models could account for the diffusional limitations, the 
network model could be more appropriate to address these issues. The one-step models such as KAS and FWO were successfully 
correlated in a single equation since both methods' difference lies in the approximation procedure for the same equation. 


Nomenclature 


A =Area [m 

Cp =Specific heat capacity [J/g.K] 

k =Thermal conductivity [W/m.K] 
L =Length [m] 

m =Mass [g] 

a =Fraction of solid materials decomposed [-] 

n =The reaction order [-] 

Wo =Initial weight of biomass [g] 

VA =weight of biomass at any run time [g] 

k =The rate constant [s+] 

ko =Frequency factor [s7}] 

E =The activation energy [k] /kmol] 
R =The universal gas constant [k] /kmol] 
YA =The reaction temperature [K] 

b =Heating rate [K/s] 

Ei =Characteristic temperature [K] 
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